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Abstract 



The theory of the 2kp and Akp instabilities in quarter-filled band organic 
conductors is revisited. The phase angles of the 2kp bond and charge density 
waves are shown to change as electron correlation is turned on, and this 
switching of the phase angle is critical for understanding the bond distortion 
patterns in the real materials. Intersite Coulomb interactions in the real 
materials must be nonzero but less than a critical value. Both intersite and 
intrasite charge density waves are destabilized in the quasi-two-dimensional 
regime for realistic parameters, thus explaining the weakening of these phases 
in the superconducting materials. 
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Conducting organic charge-transfer solids are of interest as highly doped Mott-Hubbard 
semiconductors. Nearly forty of these materials are superconducting, and evidence exists 
for strong Coulomb interactions among the fermions in these systems [1]. For a complete 
understanding of their normal state, it is essential that the spatial broken symmetries in 
the quasi-one-dimensional(quasi-one-d) nonsuperconducting conductors be precisely under- 
stood. In addition to the usual 2kp Peierls instability, many of the organic conductors exhibit 
a Akp instability (kp is the one-electron Fermi wavevector) [2]. In spite of considerable theo- 
retical work on the origin of the Akp instability, important issues involving distortions in the 
real materials and the parameter space in which the materials lie remain largely unresolved. 
This is primarily due to the limitations of calculations of susceptibilities [3-5] that do not 
measure actual distortion patterns or charge densities. We give here a considerably more 
detailed picture of the transitions in the actual experimental systems. The present picture 
is shown to have significant implications for the quasi-2-d superconductors, in which these 
transitions are absent or weakened. 

We will limit ourselves to p = 0.5, where p is the number of electrons per site. For this 
commensurate case, the bond order wave (BOW), with periodic modulation of the intersite 
distances, and the charge density wave (CDW), with modulation of the intrasite charge 
densities, are distinct. Each of these can have periodicities 2kp (period 4) and Akp (period 
2). Furthermore, each period 4 density wave can occur in two forms, corresponding to dif- 
ferent phase angles. In Table 1 we show the various period 4 and period 2 patterns. For zero 
Coulomb interactions, the ground state consists of coexisting 2k p BOW1 and CDW1 [6]. Ta- 
ble 1 also shows the experimentally observed bond distortion pattern [7] in M E M (T C N Q) 2 
below the 2k p transition. Our discussions are not limited to MEM(TCNQ) 2} rather, this 
is one system whose bond distortion pattern is precisely known. 

Existing theories [3-5,8-10] do not explain several of the mysteries in the real materials. 
Observation of signatures of electron-molecular vibration (e-mv) coupling below the Akp 
transition [11] has in the past led to the interpretation of the Akp phase as a CDW [9,10]. 
This would imply (qualitatively speaking) alternate occupancies of sites by electrons (see 
Table 1, row 1, column 3). Previous work has shown that the Akp and 2k p CDWs do not 
coexist [4]. Thus in such a case, for the 2k p transition to occur on the same chain, the 
2kp phase can only be a spin density wave (SDW) or a spin-Peierls state. We discount the 
possibility of the SDW, as long range SDW in one dimension is not possible. The spin- 
Peierls state that can accompany the Akp CDW is the BOW1 state of Table 1, in which the 
distances between the "occupied" sites are alternating. The experimental hond tetramerized 
phase, shown in the last column of Table 1, is, however, different from BOW1, and therefore 
we conclude that the Akp cannot be a CDW to begin with. The Akp phase can also be 
a BOW (Table 1, column 3, second row), and indeed in this case peaks in the wavevector 
dependent charge-transfer susceptibility are found both at 2kp and Akp within the extended 
Hubbard model [4]. What is still not clear, however, is how the experimental bond distortion 
pattern is obtained below the 2k p transition, as no combination of the Akp BOW and the 
2kp BOW1 (which dominates over the BOW2 for noninteracting electrons) can give the 
observed distortion pattern. The strong role of the e-mv coupling [11] is also diffcult to 
understand if Akp is a BOW. 

We show in the present work that these problems can be resolved if, (a) phonons are 
included explicitly in the calculations, and (b) actual bond distortion patterns and charge 
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modulations are measured. The most significant findings of our calculations are: (i) the 
phase angles of the lowest energy period 4 structures change with nonzero Coulomb inter- 
action. This switching of the phase angles is central to the detailed understanding of the 
Akp and the 2kp transitions, (ii) Exactly as the 2kp BOWf and 2kp CDWf coexist for all 
e-mv and electron-acoustic phonon (e-ph) couplings [6], the 2kp BOW2 and and the 2kp 
CDW2 also coexist, (iii) The absolute ground state, for the intersite Coulomb interaction 
less than a critical value, has a Akp BOW component, in addition to the BOW2 and CDW2 
components. The experimental lowest energy bond distortion pattern now emerges natu- 
rally as a combination of the period 2 BOW and the period 4 BOW2. (iv) Only for large 
intersite Coulomb interaction the ground state is a Akp CDW, which does not describe the 
experimental systems. 

The Hamiltonian we consider is, 

H =J2( t ~ a ( U i+l ~ U i)) (4x C * + l,<x + 4+1,<t C «t) 

id 

+ U n^n^ + V niUi+i 

i i 

K 2 K 2 

i i i 

Here c\ a creates an electron of spin a at site z, n 8(7 = c] a Ci a , = n iui t is the one-electron 
hopping integral, and U and V the on-site and nearest neighbor Coulomb repulsions. The 
Ui are the displacements of the ith molecular units from their equilibrium positions, V{ 
corresponds to an internal molecular mode, a and (3 the e-ph and e-mv coupling constants, 
and Ki and K 2 are the corresponding elastic constants. 

Exact numerical calculations were done for finite rings of 8, 12 and 16 sites for nonzero 
U and V, both with and without the restriction of constant overall ring size. We first show 
the results for the case where the overall system size is constrained. We begin with (3 = 
and initially assume the functional form for Uj that is valid in the U = V = limit, viz., 
Uj = u cos(2kpj — 9). The calculations are done for both 6 = and 6 = 7r/4, corresponding 
to the BOW1 and BOW2. The 2kp BOW is weakened by U. However, the more significant 
result is the crossover from 6 = to 6 = 7r/4, corresponding to a transition from 2kp BOW1 
to 2kp BOW2, for nonzero U. This crossover in the phase of the 2k p BOW, which we show 
is critical for a complete understanding, is shown in Fig. 1(a). 

In addition to causing a change in 6, Coulomb interactions also add a 4^-component 
to the ground state wavefunction. We repeat the calculation of the ground state energies 
with the Ansatz u } = u [r 2 cos(2kpj — 6 2 ) + r 4 (cos4:kpj — # 4 )], where r 2 + r 4 = 1. Numerical 
calculations indicate that 6 2 = 7r/4 and 64 = 0. The absolute lowest energy within correlated 
models occurs for r 4 7^ 0. Furthermore, r 4 increases monotonically with U } and for fixed U 
increases rapidly with V, provided V is less than a critical value V C (U). We discuss the issue 
of V c below. The effects of U and V on r 4 are summarized in Fig. 1(b). 

We now emphasize the following point. The experimentally observed bond distortion 
pattern (last column of Table 1) is exactly what is obtained with our Ansatz for r 4 > 0.41. 
We conclude then that the crossover in 6 and nonzero r 4 > 0.41 are both necessary to explain 
the experimental ground state bond distortion. Calculations of susceptibilities [3-5] do not 
distinguish between BOW1 and BOW2, and also cannot give the exact bond distortion 
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pattern. 



The issue of V c can be understood as follows. In the limit of U — > oo and (3 = 0, Eq. 1 
reduces to the half-filled band spinless fermion Hamiltonian, 



Transition to a CDW with equal bond lengths occurs for V > V c = 2t and small a. For 
V < V c a dimerized BOW appears for infinitesimal a, with the bond dimerization increasing 
with V [12]. Numerical simulations in the limit of a = indicate that even for finite U, a 
CDW appears only above a V C (U), which is now slightly larger than 2t [13]. The increase of 
r 4 with V within our model is then expected for V < V C (U). Thus the observed periodicity 
in the last column of Table 1 indicates that experimental materials are in the V < V c regime. 

The above results were confirmed by a second set of calculations. Instead of using the 
size constraint and the Ansatz for it 4 -, the ground state energy was minimized with respect 
to a position dependent ti using the Lanczos method in an iterative manner to find the final 
ti. These calculations also find that for U = V = the BOW1 is lower in energy, but for 
nonzero U the BOW2 is favored. 

As in the previous case, the absolute ground state is once again found to be a composite 
of the Akp BOW and the 2kp BOW2. To display the distortion patterns we compute 
the Fourier transform T(q) of the final hopping parameters ti which minimize the energy. 
Period four distortions appear as a peak in the Fourier transform at momentum 2kp } and 
the difference between the BOW1 and the BOW2 patterns is the phase of this Fourier 
component. Similarly the period 2 component of the distortion appears at momentum Akp. 
Crossover occurs from a dominant 2kp to a dominant Akp when T(2kp) = T(4&p). Our 
results are summarized in Fig. 2(a) and (b), where the 2kp curves correpond to the BOW2 
phase, which is lower in energy. As in Fig. 1(b), here also nonzero U enhances the Akp 
contribution to the wavefunction, but the Akp dominates over the 2k p only for nonzero 
moderate V . Once the Akp contribution to the wavefunction dominates over the 2kp, the 
overall bond distortion pattern is exactly the same as the experimental pattern. 

In the above we have discussed only the BOW. The experimental signature of e-mv 
coupling [11] requires also a CDW. Our numerical work finds that exactly as the BOW1 
coexists with CDW1 [6], the BOW2, and the composite ground state of BOW2 and the 
Akp BOW, coexist with the CDW2. The CDW2 is confirmed by monitoring the actual site 
charge densities. As expected, the2kp CDW2 vanishes for V > V C (U) } where the Akp CDW, 
not relevant for the experimental systems, appears. 

The above results give a complete picture of the 2k p and Akp instabilities in the real 
materials. The free energy of the experimental system is dominated by high spin states 
at high temperatures, whose electronic behavior is similar to that of the spinless fermion 
Hamiltonian of Eq. 2. For V < V c the lattice dimerizes to the Akp BOW below the Akp 
transition temperature. At still lower temperatures, the free energy is dominated by low spin 
states, and the behavior now would resemble that of the ground state of Eq.l. Dimerization 
of the dimerized lattice now takes place, and the overall bond distortion pattern resembles 
that in the last column of Table 1. Coexistence of these states with the 2k p CDW2 indicates 
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that the overall ground state is a composite of the Akp BOW, 2kp BOW2, and 2kp CDW2. 
Signatures of e-mv coupling [11] are then expected for any nonzero f3. 

Although the present work has focused on the (quasi-l-d) nonsuperconductors, our re- 
sults have strong implications for the quasi-2-d superconducting TMTSF and the BEDT- 
TTF based materials. The tendency to have a phase transitions discussed here is consid- 
erably weakened in the superconductors [1]. In the quasi-2-d regime, the Akp CDW with 
alternate occupied sites is still possible. However, our demonstration that organic conduc- 
tors have V < V c would explain the absence of this CDW transition. On the other hand, the 
absence of the BOW transition is a result of two- dimensionality. The BOW transition here 
is analogous to the spin-Peierls transition in p = 1, which is weakened in 2-d, and gives way 
to antiferromagnetism [14]. The CDW2, which we find here, is a result of the BOWs, and is 
thus also expected to vanish for sufficiently strong two-dimensional electron hopping. The 
only remaining possible transition is then to the 2k p SDW, and it has been argued [15] that 
for p away from 1 the SDW requires a minimal 2-d hopping, but then gradually weakens as 
the extent of two-dimensionality increases. This would be supported by the occurrence of 
a spin-Peierls transition in the quasi-l-d TMTT materials, the occurrence of SDW in the 
weakly two-dimensional TMTSF, its vanishing under pressure, and the absence of both spin- 
Peierle and SDW transitions in the less anisotropic BEDT-TTF. Thus all possible spatial 
broken symmetries are considerably weakened in p = 0.5 in the quasi-2-d regime for V < V c . 
Whether or not superconductivity is related to the suppression of spatial broken symmetries 
is an intriguing question. 

Some of these calculations were done on the iPSC/860 and Paragon at the San Diego 
Supercomputer Center. This work was supported by DOE grant DE-FG02-85ER-40213. 



5 



REFERENCES 



[1] See The Physics and Chemistry of Organic Superconductors, G. Saito and S. Kagoshima, 

eds. (Springer- Verlag, Berlin, 1990). 
[2] J. P. Pouget and R. Comes, in Charge Density Waves in Solids, edited by L.P. Gorkov 

and G. Gruner, (North-Holland, 1989), pp. 85-136 and references therein. 
[3] V.J. Emery in Highly Conducting One-Dimensional Solids, edited by J.T. Devreese et 

al., (Plenum, N.Y., 1979), p.247. 
[4] J.E. Hirsch and D.J. Scalapino, Phys. Rev. Lett. 50, 1168 (1983); Phys. Rev. B 27, 

7169 (1983); ibid, B 29, 5554 (1984). 
[5] S. Mazumdar, S.N. Dixit and A.N. Bloch, Phys. Rev. B 30, 4842 (1984). 
[6] K.C. Ung, S. Mazumdar and D.K. Campbell, Solid St. Commun. 85, 917 (1992). 
[7] R.J.J. Visser, S. Oostra, C. Vettier and J. Voiron, Phys. Rev. B 28, 2074 (1983). 
[8] J. Hubbard, Phys. Rev. B 17, 494 (1978). 

[9] E.M. Conwell, A.J. Epstein and M.J. Rice in Lecture Notes in Physics, Vol. 95, edited 
by S. Barisic and J.C. Cooper (Springer, Berlin, 1979), p. 204; E.M. Conwell and LA. 
Howard, Phys. Rev B31, 7835 (1985). 
[10] J. Kondo and K. Yamaji, J. Phys. Soc. Jpn. 43, 424 (1977). 

[11] R.P. McCall et al., Mol. Cryst. Liq. Cryst. 120, 59 (1985) and references therein. 
[12] P.P. Bendt, Phys. Rev B 30, 3951 (1984) and references therein. 
[13] H.Q. Lin, E. Gagliano, D.K. Campbell and J.E. Gubernatis (unpublished). 
[14] S. Mazumdar, Phys. Rev. B 36, 7190 (1987); S. Tang and J.E. Hirsch, Phys. Rev. B 
37, 9546 (1988). 

[15] S. Mazumdar, H.Q. Lin and D.K. Campbell, Synth. Metals, 41-43, 4047 (1991). 



6 



10 

9 
8 
7 
6 
5 
4 
3 
2 
1 




j i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i_ 













BOW and CDW patterns in p=0.5 






Period 4 


Period 4 


Period 2 


MEM(TCNQ) 2 






2k B0W1 

F 


2k B0W2 

F 


4k BOW 

F 


BOW 






2k CDWl 

F 


2k CDW2 

r 


4k CDW 

F 


CDW 













~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 

0123456789 10 



TABLE I. The possible bond distortions and charge modulations in a p = 0.5 chain. The double 
and dotted bonds correspond to bond lengths shorter and longer than the average (undistorted) 
bond length, while the single bonds correspond to undistorted bonds. In the last column we have 
given the bond distortion pattern in MEM(TCNQ)2 below the 2kp transition. Here the double 
dotted bond is not as long as the single dotted bond. In case of the CDWs, the lengths of the 
vertical bars on the sites correspond to the charge densities. 
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(a) U (b) U 

FIG. 1. The gain in total energy vs. U for 2k p bond distortions with two different phase angles 
for N = 8 and 12. (a) (N is the length of the system.) The solid line in each case corresponds to 
the BOWI and the dashed line to the BOW2. The arrows indicate the crossing of the BOW1 and 
BOW2 energies, (b) The relative weight r± of the Akp component of the ground state wavefunction 
as a function of U and V. The dashed and solid curves correspond to N = 8 and 12 respectively. For 
r4 > 0.41 the bond distortion pattern is that observed in MEM(TCNQ) below the 2kp transition. 
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(a) U (b) V 



FIG. 2. The 2kp (solid lines) and 4kp (dashed lines) components of the hopping parameters 
at V = and A = a 2 / K\ = 0.45. (a) Results are shown for N = 8 (circles) and N = 16 
(diamonds) with periodic boundary conditions (circles), and for N = 12 with antiperiodic boundary 
conditions (squares). In all cases we see that the 4kp component increases with U until at very 
large U the distortion vanishes completely. Note that, U alone is insufficient to drive the 4k p 
component larger than the 2kp component, which is the condition for the pattern to change from 
"medium-short-medium-long" to "short-medium-short-long", (b): Same as (a) for fixed U = 8 
and varying V. Note that V enhances the 4kp component, which can now be larger than the 2kp 
component. The results for L = 16 at large V may be an artifact of the small system size. At the 
end of this curve, &tV = 6, the L = 16 distortion is purely 4kp, or period two. However, the wave 
function at this point has period of four, changing sign under translation by two lattice sites. For 
the undistorted lattice, this corresponds to a ground state with nonzero momentum. 
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